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In this paper, the operation of totally asymmetric simple exclusion process with 
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. mathematical analysis and numerical simulations, we have found that, according to 
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the method chosen by the particle at the bifurcation, the model can be separated 
into two different situations which lead to different results. The results obtained in 
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I. INTRODUCTION 



y asymmetric simple exclusion process (TASEP) has been 

m 

. It is one special type of ASEPs which represent one of 



One-dimensional total 
studied for many years 
the basic models studying non-equilibrium behavior of particle transport along one- 
dimensional lattices. ASEP was first introduced in [Sj for the explanation of ribosome 
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inside the cell of creature. Now, ASEP has also been used to simulate a lot of physical 
p.oce.ses .c.„di,, surface d,ffus,o. ^J, traffic .ode. g and .oleo.la. ,noto. fi.Q, 
etc. 

TASEP, as a special case of ASEP, has two identical characters. One is that the 
model is discrete, which means that it is studied in finite time intervals. The other 
is that the particles in the lattices can only move in one direction. The analytical 
solutions under open boundary conditions had been obtained in 

Under open boundary conditions, the solutions yield phase diagrams with three 
phases jsl, O]. At small values of injection rates a < 0.5 and a < [3, the system is 
found in a low-density entry-limited phase where 

Pi = a pl = J = a{l-a) Phuik = a (1) 

where pi, pi and p^juik are the densities at the entrance, exit and the bulk of the lattice 
far away from the boundaries, respectively. J denotes the flux. 

At small values of extraction rates /3 < 0.5 and (3 < a, the system is in a 
high-density exit-limited phase with 

Pi = l- ^^^~ PL = l-(3 J = (3{l-P) p,uik = l-(3 (2) 

a 

At large values of the injection (a ^ 0.5) and extraction {(3 ^ 0.5) rates the system 
is in a maximal current phase with 



Pi = l-^ PL = ^ J=\ P^ik = \ (3) 
A large number of varieties of extensions of TASEPs have been investigated, such as 



TASEP with 
speed TASEP 



lierarchical long-range connections {lo| . two lane situations 



111, two 



12| . the effect of defect locations[13], particle-dependent hopping rates 
iJl, and so forth. 

n ^ 

Recently, J. Brankov pj| and E. Pronina [16] studied an ASEP with two chains 
in the middle of the filament. They supposed that a particle chooses to move into 
these two chains with equal probability, 0.5 and these two chains have the same 
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length. Then in 2007, Yao-Ming Yuan and Rui Jiang [l^] investigated a TASEP with 
a shortcut in the middle. They set a possibility q for a particle to jump through the 
shortcut when it faces the bifurcation, and they set the length of the shortcut to be 
zero. Unfortunately, they did not state the situation clearly. 

In this paper, we will investigate TASEP with one or two shortcuts respectively. 
The length of the shortcut is also assumed to be zero. For example, the filament on 
which the motor moves may be twisted as figure 1(a), a motor may have a chance to 
jump directly from site ki to A;2, as shown in figure 1(b). 

To the basic model, i.e. there is only one shortcut along the filament. It is natural 
for us to divide the whole filament into three segments. As shown in figure 1(b), a 
molecular motor will face a choice of whether to jump through the shortcut or to 
move ordinarily through segment 2, when it reaches site ki. An important problem 
is that the particle at site ki may have to wait to go through the shortcut if there 
is also another particle occupying site k2 — 1- Which of them to go first should be 
determined clearly. 

After the basic model investigated, we would like to do some research work on our 
advanced model 1. In advanced model 1, there are two shortcuts which begin and 
end at different sites. These two shortcuts stay respectively along the filament, that 
is to say, shortcut 2 begins after shortcut 1 ends, as shown in figure 2. 

Also, considering the difference of places where the shortcuts stay, we have ad- 
vanced model 2, as shown in figure 3. In advanced model 2 there are also two short- 
cuts. But they begin at the same site and end at different, which means when a motor 
moves into site ki, it will choose whether to pass through shortcut 1 or shortcut 2 or 
to go ordinarily through segment 2. 

In the following, we firstly give the detailed description of the different models, 
then we will discuss their phase situation of the corresponding segments theoretically, 
which is followed by the numerical simulations. 
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(a) 




ki+l 



(b) 



FIG. 1: basic model with one shortcut: (a) a road with one shortcut at the bulk; (b) lattices with 
one shortcut. 




(a) 




(b) 



FIG. 2: advanced model with two shortcuts at different beginnings and endings: (a)a road with two 
shortcuts at different beginnings and endings; (b)lattices with two shortcuts at different beginnings 
and endings. 
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(a) 




(b) 



FIG. 3: advanced model with two shortcuts at the same beginning and different endings: (a)a road 
with two shortcuts at the same beginning and different endings; (b)lattices with two shortcuts at 
the same beginning and different endings. 

II. MODELS 
A. Basic model 



First, we divide the lattices into three segments, as shown in figure 1(b). Segment 
1,2,3 start from site 1, A;i + 1, ^2 and end in site ki,k2 — 1, L, respectively. 

In this model, all the particles traveling in the lattices are identical to each other. 
At time t = 0, we may suppose that the lattices are empty. As time goes, there will be 
a particle in site 1. This particle will pass continuously through the whole filament. 

Given a site i(l ^ z ^ L), in an infinitesimal time interval dt, ii i = 1, a particle 
is inserted with probability adt, provided the site is empty, if i = L and there is 
a particle in it, then the particle in site L is extracted with probability /3dt, and if 
1 ^ i ^ L{i ^ ki, k2 — I) and site i is occupied, the particle in site i will move into 
site i + 1 with probability dt, providing site i + 1 is empty. However, if a particle is in 
site ki or k2 — 1, we have to clarify the way it moves. Although Yao-Ming Yuan,Rui 
Jiang have discussed the method of choice at site ki, the situation at site k2 — I 
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has not been studied clearly. It is possible that when both site A;i and site k2 — l have 
particles, it may produce a collision in site k2- 

1. Situation 1 

If site ki has a particle, then the particle will choose the way as following: 

• if site ki-\-l and site ^2 are both occupied, then the particle in site k\ does not 
move. 

• if site ^2 is occupied but site fci + 1 is empty, then the particle in site ki will 
move into site /ci + 1 with probability dt. 

• if site /ci + 1 is occupied but site k2 is empty then: 

•k if site A;2 — 1 is empty, the particle in site k\ will move into site k2 with probability 
qdt. 

• if site k2 — 1 is occupied , the particle in site ki will move into site k2 with 
probability qpdt. 

• if both site k\-\- 1 and site ^2 are empty, then the particle in site ki will move 
into site ki-\-l with probability (1 — q)dt, and: 

■k if site A;2 — 1 is empty, the particle in site ki will move into site /c2 with probability 
qdt. 

• if site A;2 — 1 is occupied , the particle in site ^1 will move into site k2 with 
probability qpdt. 

If site k2 — 1 has a particle, and site k2 is empty, then the particle will choose the 
way as following: 

• if site ki is empty, then the particle in site A;2 — 1 will move into site k2 with 
probability dt. 

• if site ki is occupied, then the particle in site ^2 — 1 will move into site k2 with 
probability (1 — pq)dt. 

The character of this model is that the particle in site ki first choose whether to 
move through the shortcut and once it has decided to jump through the shortcut, it 
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will face the problem of whether site A;2 — 1 is occupied. If site A;2 — 1 is occupied, then 
the particle which has chosen the shortcut can only move into site k2 with probability 
pdt in order to prevent collision. This means that the particle takes a risk to jump 
the shortcut because once it has chosen the shortcut it can not go back into segment 
2, even if site /ci + 1 is empty. 

2. Situation 2 

If site ki has a particle, then the particle will choose the way as following: 

• if site ki + 1 and site k2 are both occupied, then the particle in site ki does not 
move. 

• if site k2 is occupied but site A;i + 1 is empty, then the particle in site ki will 
move into site ki + I with probability dt. 

• if site fci + 1 is occupied but site k2 is empty then: 

ic if site A;2 — 1 is empty, the particle in site ki will move into site /c2 with probability 
qdt. 

■*r if site A;2 — 1 is occupied , the particle in site ki will move into site k2 with 
probability qpdt. 

• if both site ki + 1 and site /c2 are empty, then: 

ic if site A;2 — 1 is empty, the particle in site ki will move into site ^2 with probability 
qdt, and it will move into site ki + 1 with probability (1 — q)dt. 

■*r if site A;2 — 1 is occupied , the particle in site ki will move into site k2 with 
probability qpdt, and it will move into site ki + 1 with probability (1 — qp)dt. 

If site k2 — 1 has a particle, and site k2 is empty, then the particle will choose the 
way as following: 

• if site ki is empty, then the particle in site ^2 — 1 will move into site ^2 with 
probability dt. 

• if site ki is occupied, then the particle in site k2 — 1 will move into site k2 with 
probability (1 — pq)dt. 

The character of this model is that once the particle has chosen the shortcut but 
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see that the particle in site k2 — 1 has priority to move into site k2, it can change its 
mind and return to site ki + 1, providing site /ci + 1 is empty. This model is closer 
to reality because actually bus drivers can change their mind in time when they are 
told that the road ahead is to some extent blocked. It is not reasonable for drivers to 
take a shortcut which is actually difficult to go through. 

B. Advanced model 1 

Now we divide the lattices into five segments, as shown in figure 2(b), segment 
1,2,3,4,5 starts in site l,ki + I,k2,k3 + l,k4 and ends in site ki,k2 — Ijk^jk^ — 
1, L, respectively. 

In this model, all the particles traveling in the lattices are identical to each other. 
At time t = 0, we may suppose that the lattices are empty. As time goes, there will 
be a particle goes into site 1. This particle will pass continuously through the whole 
filament. 

Given site i(l ^ i ^ L), in an infinitesimal time interval dt, ii i — 1, a particle 
is inserted with probability adt, provided the site is empty, ii i — L and there is 
a particle in it, then the particle in site L is extracted with probability f3dt, and if 
1 ^ i ^ L{i 7^ ki, k2 — l, ks, k4^ — l) and site i is occupied, the particle in site i will move 
into site i + 1 with probability dt, providing site i + 1 is empty. However, if a particle 
is in site ki or A;2 — 1 or ks or k4 — 1, we have to clarify the way it moves. We will 
construct advanced model 1 using similar method according to basic model(situation 

!)• 

If site ki has a particle, then the particle will choose the way as following: 

• if site ki + 1 and site ^2 are both occupied, then the particle in site ki does not 
move. 

• if site ^2 is occupied but site /ci + 1 is empty, then the particle in site ki will 
move into site ki + 1 with probability dt. 

• if site /ci + 1 is occupied but site k2 is empty then: 

~k if site A;2 — 1 is empty, the particle in site ki will move into site k2 with probability 



9 



Qidt. 

ic if site A;2 — 1 is occupied , the particle in site ki will move into site k2 with 
probability qipidt. 

• if both site ki + 1 and site ^2 are empty, then the particle in site ki will move 
into site ki + 1 with probability {1 — qi)dt, and: 

ic if site A;2 — 1 is empty, the particle in site ki will move into site k2 with probability 
Qidt. 

•k if site ^2 — 1 is occupied , the particle in site ki will move into site k2 with 
probability qipidt. 

If site k2 — 1 has a particle, and site ^2 is empty, then the particle will choose the 
way as following: 

• if site ki is empty, then the particle in site k2 — 1 will move into site k2 with 
probability dt. 

• if site ki is occupied, then the particle in site k2 — 1 will move into site k2 with 
probability (1 — piqi)dt. 

If site k^ has a particle, then the particle will choose the way as following: 

• if site /cs + 1 and site /C4 are both occupied, then the particle in site k^ does not 
move. 

• if site ki is occupied but site A;3 + 1 is empty, then the particle in site k^ will 
move into site ^3 + 1 with probability dt. 

• if site /cs + 1 is occupied but site A;4 is empty then: 

-k if site A;4 — 1 is empty, the particle in site k^ will move into site /C4 with probability 
q2dt. 

ic if site A;4 — 1 is occupied , the particle in site ks will move into site ki with 
probability q2P2dt. 

• if both site ^3 + 1 and site ki are empty, then the particle in site ks will move 
into site A;3 + 1 with probability (1 — q2)dt, and: 

•k if site A;4 — 1 is empty, the particle in site k^ will move into site k^ with probability 
q2dt. 

• if site A;4 — 1 is occupied , the particle in site kz will move into site ki with 
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probability q2P2dt. 

If site ki — 1 has a particle, and site is empty, then the particle will choose the 
way as following: 

• if site /c3 is empty, then the particle in site — 1 will move into site k^ with 
probability dt. 

• if site ks is occupied, then the particle in site k4 — 1 will move into site A;4 with 
probability (1 — P2Q2)dt. 

C. Advanced model 2 

Also we can divide the lattices into four segments, as shown in figure 3(b), segment 
1,2,3,4 starts in site 1, ki + 1, k2, ks and ends in site ki, k2 — l,ks — 1, -L,respectively. 

In this model, all the particles traveling in the lattices are identical to each other. 
At time t = 0, wc may suppose that the lattices are empty. As time goes, there will 
be particles goes into site 1. These particles will pass continuously through the whole 
filament. 

Given site i(l ^ i ^ L) , in an infinitesimal time interval dt, if i — 1, a, particle 
is inserted with probability adt, provided the site is empty, if i = L and there is 
a particle in it, then the particle in site L is extracted with probability pdt, and if 
1 ^ i ^ L{i ^ fci, ^2 — 1, ^3 — 1) and site i is occupied, the particle in site i will move 
into site i+1 with probability dt, providing site i + 1 is empty. However, if a particle 
is in site A;i or ^2 — 1 or — 1, we have to clarify the way it moves. We will construct 
advanced model 2 using similar method according to basic model(situation 1). 

If site ki has a particle, then the particle will choose the way as following: 
• If sites ^1 + 1,A;2 and ks are all occupied, the the particle in site ki does not move. 

• If site A;i + 1 is empty but sites k2 and k^ are occupied, the particle in site ki 
will move into site ki + I with probability dt. 

• If site k2 is empty but sites ki + 1 and k^ are occupied, then: 

ic If site ^2 — 1 is empty, the particle in site ki will move into site k2 with probability 
qidt. 
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-k If site ^2 — 1 is occupied, the particle in site ki will move into site k2 with 
probability qipidt. 

• If site fcs is empty but sites ki + 1 and k2 are occupied, then: 

~k If site A;3 — 1 is empty, the particle in site ki will move into site k^ with probability 
q2dt. 

~k If site — 1 is occupied, the particle in site ki will move into site ks with 
probability q2P2dt. 

• If site ki + 1 and site k2 are empty but site k^ is occupied, then the particle will 
move into site ki + 1 with probability {1 — qi)dt, and: 

• If site ^2 — 1 is empty, the particle in site ki will move into site k2 with probability 
qidt. 

• If site A;2 — 1 is occupied, the particle in site ki will move into site k2 with 
probability qipidt. 

• If site ki + 1 and site k^ are empty but site k2 is occupied, then the particle will 
move into site ki + 1 with probability {1 — q2)dt, and: 

-k If site — 1 is empty, the particle in site ki will move into site ks with probability 
q2dt. 

If site /cs — 1 is occupied, the particle in site ki will move into site k^ with 
probability q2P2dt. 

• If site k2 and site k^ are empty but site A;i + 1 is occupied, then: 

~k If site A;2 — 1 is empty, the particle in site ki will move into site k2 with probability 
qidt. 

~k If site A;2 — 1 is occupied, the particle in site ki will move into site k2 with 
probability qipidt. 

• If site — 1 is empty, the particle in site ki will move into site ks with probability 
q2dt. 

•k If site /cs — 1 is occupied, the particle in site ki will move into site k^ with 
probability q2P2dt. 

• If sites ki + 1, k2 and ks are all empty, then the particle in site ki will move into 
site ki + 1 with probability {1 — qi — q2)dt, and: 



12 



•k If site — 1 is empty, the particle in site ki will move into site k2 with probability 
qidt. 

If site /c2 — 1 is occupied, the particle in site ki will move into site /c2 with 
probability qiPidt. 

~k If site k^ — l is empty, the particle in site ki will move into site ks with probability 
q2dt. 

• If site fcs — 1 is occupied, the particle in site ki will move into site ks with 
probability q2P2dt. 

If site A;2 — 1 has a particle, and site k2 is empty, then the particle will choose the 
way as following: 

• if site ki is empty, then the particle in site k2 — 1 will move into site k2 with 
probability dt. 

• if site ki is occupied, then the particle in site k2 — I will move into site k2 with 
probability (1 — piqi)dt. 

If site ks — 1 has a particle, and site ks is empty, then the particle will choose the 
way as following: 

• if site ki is empty, then the particle in site /ca — 1 will move into site k^ with 
probability dt. 

• if site ki is occupied, then the particle in site ks — 1 will move into site ^3 with 
probability (1 — p2q2)dt. 

III. ANALYSIS OF MODELS 

A. Theoretical Analysis of the Basic Model 

In this section, we will give a theoretical analysis of the phase situation of three 
segments of our basic model. We choose situation 1 of our basic model described 
above as the base of our analysis. 

According to figure 1(b), we are able to express the insertion and extraction rates 
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using densities of certain sites. These rates are the following: 

tte//! = a 

Peffl = (1 - Pfel+l)Pfe2 + (1 - Pfe2)(l - Pfel+l)(l -Q)+ ?(1 - Pfe2)[(l - Pk2-l) + PPk2-l] 

aeff2 = PkAPk2 + (1 - Pfe2)(l - ?)] 

/3e//2 = (1 - Pfc2)[(l - Pfci) + ?(1 -p)Pfci + (1 - 5)pfei] 
ae//3 = Pfei(l - Pk2-l)Q + Pfe2-l(l - Pfel) + PfeiPfe2-l 
/3e//3 = /3 

</sC = QPfcl(l - Pfc2)[l - Pfc2-1 +PPfc2-l] 

(4) 

Certainly, the flux should be conserved: 

Jl — J2'^ Jsc — (5) 

Firstly, we are going to prove that segment 1 and segment 3 should be in the same 
phase. 

If segment 1 is in maximum-current phase, then we have Ji — \ according to (3). 

According to (5), wc have J3 = |. So segment 3 is also in maximum-current phase. 
Secondly, we claim that the following four cases are impossible 

(L, L, H), (L, H), {H, L, L), {H, H, L) 

To prove that, we firstly focus on the first two situations. If the phase condition is 
(L, *, H), due to (1), we have: 

Pi^a pL^l-P, {a,/3<0.5) (6) 

Because Ji = J3, we obtain: 

a(l - a) = (3(1 - 13) (7) 

When a ^ (3, from (7) we have a + (3 = 1, which contradicts a,P < 0.5. Thus the 
phase situations (L, *,H) are impossible. 

If the phase condition is {H,L,L), it is difficult for us to analysis the original 
model with parameter p. So we set p = 1, which means that the particle in site ki 
has priority to move into site k2 when site A;2 — 1 is occupied. 
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Due to (2), we have: 

/?e//l = 1 - Pifci = 1 - Pfei+1 + gPfei+l(l - Pfca) (8) 

from which we can obtain: 

Pfei =Pfci+i(l-g(l-Pfc2)) (9) 
Let 1 — q{l — pfcj) = ^) then (9) turns into: 

Pfei = Apk^+i (10) 

Due to (1), we have: 

ae//2 = Pfci+l = /3fei[l-g(l-pfe2)] (11) 

The equation above can be turned into: 

Pfei+i = ^Pfei (12) 

From (10) and (12), we reach: 

Pfei = ^Vfei (13) 

If = 0, then Ji = Pa:i(1 — PfcJ = does not agree with reahty. So could not 
be 0. Thus A^l. For all g € [0, 1], 

Pfc2 = 1 (14) 
On the other hand, due to (5), we have: 

Pfcl(l-Pfel) =Pfc2(l-Pfc2) (15) 

Because pk^ = 1 — Peffi > 0.5, pk^ < 0.5, we know that pk^ ^ pk^- Together with (15), 
we have p^^ + p^^ = 1. Together with (14), we get p^^ = 0. It is impossible. Thus the 
phase situation can't be {H, L, L). 

If the phase condition is {H, H, L), also we can set = 1 for our convenience to 
prove. Because segment 2 is in high-density phase, we have /3e//2 = 1 — Pfe2-i- Because 
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segment 3 is in low-density phase, we have Q;e//3 = Pk2- According to (4), we have 
the cases: 

1 - Pfc2-1 = 1 - Pfe2 - 5Pfcl(l - Pfe) 
Pfe2 = ?Pfcl(l - Pfc2-l) + Pfc2-1 

Let gpfej = then the cases become: 



(16) 



(17) 



Pfe2-1 Pfe2 + -^(1 - Pfc2) 
Pfc2 = Pfc2-1 + B{1 - Pk2-l) 

From (17) we can obtain: 

{2B - B^)pk, = {2B - B^) (18) 

Because J3 = ^^2(1 ~ Pfc2)) we know that 7^ 1, together with (18), we will have 
= 0, which leads to pk^ = 0, which does not agree with reality. Thus {H, H, L) has 
been excluded. 

Now we have proved that segment 1 and segment 3 should be in the same phase. 
In the following, we will discuss which phase may segment 2 be in. 

If the three segments are in (L, *, L), from the proof which has demonstrated that 
{H, H, L) is impossible, we know that (L, H, L) can not be true either. Thus we can 
definitely say that for (L, *, L), only (L, L, L) can be right. 

If the three segments are in (if, *, H), let us see whether [H, L, H) is true. If the 
three segments are in (if, L, ii), then we have: 

I3effl = 1 - pfei = 1 - + 9Pfci+l(l - Pfc2) 

(19) 



a 



;//2 = Pfei[l - g(l - Pfea)] = Pfei+i 



Let 1 — g(l — pfcj) = A, then from (19) we can get p^^ = Ap^^, then A = 1, then 

Pfe2 = 1- 

Because of (5), we have: 

Pfei(i - Pfei) = Pfei+i(i - Pfei+i) + gpfei(i - Pfc2) = - /5) (20) 



16 



If we put pfe2 = 1 into (20), we will get: 

- Pfei) = Pfei+i(l - Pfei+i) (21) 

Because pk^ ^ Piki+i, we have pfe^ + p^i+i = 1. From p].^ — \ — /3, we know that 
Pifci+i = /3, then J2 — P{1 — P) — J3, which does not agree to reahty. So (if, L, if) is 
impossible. 

Thus we have excluded {H, L, H). We now know that for (if, *, ii), only (if, if, if) 
is possible. 

As for (M, *, M), it is a tough task for us to analyze the situation of segment 2 
using equations. But under assistance of numerical simulations, it will be possible for 
us to find which phase is true for segment 2. 

For situation 2 of model A, we can also write the insertion rates and extraction 
rates at special sites: 

tte//! = « 

/3e//l = (1 - Pfci+l) + gPfei+l(l - Pfc2)[(l - Pfe2-l) +PPfc2-l] 

aeff2 = PkApk2 + (1 - Pfc2)(l - - Pfc2-l) +PPfe2-l))] , . 

[11) 

Peff2 = (1 - Pfc2)[(l - Pfei) + - P)pki + (1 - 9)Pfci] 
«e//3 = Pfei(l - Pfe2-l)? + Pfc2-l(l - Pfci) + PfclPfc2-l 
/3e//3 = /? 

The analysis of situation 2 of basic model is the same of the analysis of situation 1, 
providing p = 1, so we do not state it repeatedly. 

From the analysis above, we can conclude that for the situation p — 1, the three 
segments should be among the following situations: (L, L, L), (if, H, if), (M, *, M). 
We will do further investigation using numerical simulations in the following sections. 

B. Analysis of advanced model 1 

In this section, we will give the phase situation of the five segments of advanced 
model 1. According to figure 2(b), we are able to express the insertion and extraction 
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rates using densities of certain sites. These rates are the following: 

aeffi = a 

Peffl = (1 - Pfel+l)Pfe2 + (1 - Pfc2)(l - Pfel+l)(l - Ql) + - Pfe2)[(l - Pk2-l)+PlPk2 

ae//2 = PfeJPfca + (1 - Pfc2)(l - 

Peff2 = (1 - Pfc2)[(l - Pfci) + - />l)Pfei + (1 - 9l)Pfci] 

«e//3 = Pfci(l - Pfc2-l)?l + Pfe2-l(l - Pfei) + PfeiPfe2-l 

/?e//3 = (1 - Pfc3 + l)Pfc4 + (1 - Pfc4)(l - Pfc3+l)(l - + ^2(1 - Pfc4)[(l - P/C4-1) +P2pfc4 

«e//4 = Pfc3[Pfc4 + (1 - Pfe4)(l - ^2)] 

Peff4 = (1 - Pfe4)[(l - Pfe4) + ^2(1 - P2)Pfe3 + (1 - g2)Pfe3] 
aeff5 = Pfc3(l - Pfe4-l)92 + Pfe4-l(l - Pfes) + Pfe3Pfe4-l 

/^e//5 

-'sci = ?lPfel(l - Pfc2)[l - Pfe2-1 +PlPfe2-l] 
JsC2 = ?2Pfc3(l - Pfc4)[l - Pfe4-1 +P2Pfc4-l] 

(23) 

Also, the flux should satisfies: 

Jl — J2 + JsCl — J^ — Ji-\- JsC2 — J5 (24) 

In the following, we will prove that the segment 1,3 and 5 are in the same phase. 
Conclusions obtained in the proof of basic model will be used in the following steps. 
For the sake of simplicity, wc always set pi = p2 = 1. 

If one of the segments 1, 3 or 5 is in maximum-current phase, then the other two 
will also be in maximum-current phase, since Ji — J3 — J5 — \- Due to (6) and (7), 
we also can find that the situation is {L, H) is impossible. 

In the following, we discuss the situation {H, L). If segment 3 is in low- 

density phase, then segment 1,2 and 3 will be in [H, *, L). According to the proof of 
the basic model, we know that this situation is impossible. Because we can regard 
segment 1,2 and 3 as basic model, only by changing /3 to be Peffs- This change 
does not influence the proof since the rate expression of Peffs is independent of the 
densities of sites in the first three segments. 
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Similarly, segment 3 can not be in high-density phase since segment 3,4 and 5 will 
be in {H, *, L). So now we have excluded (L, *, *, *, H) and (i?, *, *, *, L),. 

On the other hand, if the phase situation is (L, *, *, *, L), then we know immedi- 
ately that segment 3 can not be in high-density phase. So segment 1,3 and 5 should 
all be in low-density phase. If the situation is {H, *,*,*,H), then we know immedi- 
ately that segment 3 can not be in low-density phase. So segment 1,3 and 5 should 
all be in high-density phase. 

Since we have proved that (L, H, L) and {H, L, H) does not exist in basic model, 
so we can conclude that all of the five segments should all be in low-density phase, 
providing segment 1 is in low-density phase, and all be in high-density phase, when 
segment 1 is in high-density phase. 

As for the situation (M, *, M, *, M), it is difficult for us to analyze it theoretically. 
So we will use numerical simulations to see whether segment 2 and 4 are in high or 
low-density phase. The simulations will be presented in the following section. 

C. Analysis of advanced model 2 

In this section, we are going to demonstrate the phase situation of four segments 
of advanced model 2. According to figure 3(b), we are able to express the insertion 
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and extraction rates using densities of certain sites. These rates are stated as follows: 

Q!e//i = a 

Peffl = (1 - Pfel+l)Pfe2Pfe3 + (1 - Pfe2)Pfel+lPfe3?l[(l - Pfe2-l) +PlPfe2-l] 

+ (1 - Pfc3)Pfci+lPfc292[(l - Pfca-l) +P2Pk3-l] 

+ (1 - + - Pfc2)Pfc3[l -Ql+ - Pfe2-l) +PlPfc2-l]] 

+ (1 - Pfci+l)(l - Pfc3)Pfc2[l - ?2 + g2[(l - Pfc3-l) +P2Pfe3-l]] 

+ Pfcl+l(l - Pfe2)(l - Pfe3)[9l[(l - Pfc2_l) +PlPfc2_l] + 52[(1 - Pfc3-l) +P2Pfe3-l]] 
+ (1 - pfei+l)(l - pfe2)(l - Pfe3)[l -qi-q2+ - Pfc2-l) +PlPfe2-l] 

+ g2[(l -Pfea-l) +P2Pfc3-l]] 
ae//2 = PkAPk2Pk3 + (1 - Pk2)Pk3i'^ - Ql) + (1 - Pfc3)Pfc2(l - ?2) 
+ (1 - Pfe2)(l - Pfc3)(l - 5l - Q2)] 

Peff2 = (1 - Pfe2)[(l - Pfel) + - Pl)Pfci + (1 - 5l)Pfei] 

ae//3 = Pfci(l - Pk2-l)qi + Pfe2-l(l - Pfei) + PklPk2-l 

Peff3 = (1 - Pfc3)[(l - Pfci) + ^2(1 - P2)pfei + (1 - Q2)Pfci] 
ae//4 = Pfei(l - Pk3-l)Q2 + Pfe3-l(l - Pfei) + PkiPk3-l 
PeffA = /3 

•/sci = giPfcl(l - Pfe2)[l - Pfc2-1 +PlPfe2-l] 
JsC2 = ?2Pfcl(l - Pfc3)[l - Pfe3-1 +P2Pfe3-l] 

(25) 

Also we have the equation of flux: 

Jl = J4 = J3 + Jsc2 ^^g-^ 

Js — J2 + Jsci 

Our conclusion is that the segments 1 and 4 are in the same phase. In order to 
simplify our proof, we always set pi = p2 = I. 

Firstly, if segment 1 or segment 4 is in maximum-current phase, then the other 
should also be in maximum-current phase due to (26). 

If the phase condition is (L, *, *, H), then due to a{l — a) — /3(1 — /?), we know 
that when a P we have a + P = 1, which is contradictory to a < 0.5 and /3 < 0.5. 
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Now we discuss the situation {H, *, *, L). If the segment 3 is in low-density phase, 
then we have: 

Jl = Piti(l - Pki) ^ J3 + Jsc2 = Pik3(l - Pka) = Ji (27) 

and 

J3 ^ Pk,{l - Pk^) (28) 
If segment 2 is in low-density phase, then 

aeff2 = Pkr = Pfci+l(l - Ql(l - Pk2) - 52(1 - Pks)) (29) 

Let 1 - qi{l - pk^) - ^2(1 - pka) = A we have: 

Pki = Apk^+i (30) 
Since segment 1 is in high-density phase, we have: 

Peffl = 1 - Pfei = 1 - Pfei+1 + giPfci+l(l - Pfea) + ?2Pfei+l(l - Pks) (31) 

So we have: 

Pfei+i = ^Pfei (32) 

From (30) and (32), we know that: 

Pki = ^Vfei (33) 

From (33) we have A = 1 because p^^ can not be 0. So pk2 = Pks = ^ for all 
?2 G (0, 1], which implies J4 = 0. So segment 2 can not be in low-density phase. If 
segment 2 is in high-density phase, then we have: 

Peff2 = 1 - Pfc2-1 = 1 - Pfe2 - ?lPfei (1 - Pfe) ^2^^ 
«e//3 = Pk2 = ?lPfci(l - Pfe2-l) + Pfc2-1 

From (34) we know that p^^ = or p^^ — 0, which will lead to Ji or J3 = that 
does not make sense when qi ^ 0. So we have excluded {H, *, L, L). 
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If the segment 3 is in high-density phase, then 

/?e//3 = 1 - Pkz-i Pk3- Q2pki (1 - Pks) 

(35) 

Olef4 = Pks = 92Pfei(l - Pfe3-l) + Pfe3-1 

Prom (35) we know that pk^ = or = 0, which will lead to Ji = or J4 = that 
does not make sense, (i?, *, H, L) is also impossible. 

Now we have proved that segment 1 and segment 4 are in the same phase. We 
now discuss the situation (L, *, *, L). Prom the proof above we know that segment 
3 can not be in high-density phase, either can not segment 2. Thus the phases can 
only be (L, L, L, L). 

While considering (ff, *, *, if), we should first assume that segment 3 be in low- 
density phase. Prom the proof above we will know that segment 2 can not be either 
in low or in high-density phase. And if we assume segment 3 to be in high-density 
phase, we know that segment 2 can not be in low-density phase. Thus the phases can 
only be {H,H,H,H). 

To the situation (M, *, *, M), it is rather difficult for us to analyze it theoretically. 
So we also use numerical simulations to get the results. 

IV. NUMERICAL SIMULATIONS OF THE DIFFERENT MODELS 

A. Simulations of basic models 

In this section, we use numerical methods to simulate basic models. Pirst we 
demonstrate the results of situation 1 and situation 2. 

According to analysis in last section, we set insertion rate a and extraction rate (3 
to be 0.3 and 0.8 for (L, L, L), 0.8 and 0.3 for {H, H, H), 0.8 and 0.8 for (M, *, M). 
The results of the simulations are plotted in figure 4. 

Prom figure 4, we can see that the results of numerical solutions are in agreement 
with the results of equation analysis. When a = 0.3 and P = 0.8, all the three 
segments are in low-density phase. When a = 0.8 and /3 = 0.3, all the three segments 
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FIG. 4: Numerical simulations of basic model (situation 1) corresponding to different phases: (a)a = 
0.3, /3 = 0.8; (b)a = 0.8, /? = 0.3; (c)a = /3 = 0.8, initial densities of all sites are high; (d)a = (3 = 0.8, 
initial densities of all sites are low. 
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are in high-density phase. When a — (3 — 0.8, segment 1 and segment 3 are in 
maximum-current phase. 

From figure 4(a), we can find that the density in segment 2 decreases when q 
increases, providing p is constant. For the same q, when p increases, the density in 
segment 2 decreases. 

From figure 4(b), we see that the density in segment 2 increases when q increases, 
providing p is constant. For the same q, when p increases, the density in segment 2 
increases. This is because when p is relatively large, more particles are able to jump 
through the shortcut to reach site k2, thus leaving particles in segment 2 in a serious 
traffic jam. 

From figure 4(c) and (d), we see that when segment 1 and segment 3 are both 
in maximum-current phase, segment 2 can be either in high-density phase or low- 
density phase, depending on the initial density rate of each lattice. If the average 
initial density rate of the lattices is high, then segment 2 will be in high-density phase 
when the fiux becomes stable, if the average initial density rate of the lattices is low, 
then segment 2 will be in low-density phase when the fiux becomes stable. 

Now wc use numerical methods to simulate situation 2 of basic model. The in- 
sertion and extraction rates are set as the same as in the simulation of situation 1. 

From figure 5 we are able to see that the results of situation 2 of basic model 
seem to be similar to the results of situation 1. But after careful inspection, we may 
find that there arc several differences between the results of the two situations, which 
manifest the distinct characters of these two situations. 

In order to do further investigation, we use figures which can show the differences 
between two situations. We calculate the values which are obtained by subtracting 
the density of each lattice of situation 2 from the density of each lattice of situation 
1. When we do these subtractions, we keep parameters q and p in constant. The 
results are shown in figure 6. In figure 6(a) and 6(b), we can see that the differences 
between the density of two situations are within 3 x 10"^^, which means that when 
p — 1, there is nearly no difference between the two situations. That is because when 
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FIG. 5: numerical simulations of basic model(situation 2) corresponding to different phases:(a)a = 
0.3, P = 0.8; (b)a = 0.8, /? = 0.3; (c)a = /3 = 0.8, initial densities of all sites are high; (d)a = /3 = 0.8, 
initial densities of all sites are low. 
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FIG. 6: The density of each lattice of situation 1 of basic model minus the density of each lattice of 
situation 2 of basic model with insertion rate a — 0.3 and extraction rate /3 = 0.8: (a) when q — 0.5 
and p = 1] (b) when q = 0.9 and p = I; (c) when q = 0.1 and p = 0; (d) when q = 0.9 and p = 0. 
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p — 1, the particle at site ki{x — | in the figure) has priority to get to site k2{x — | 
in the figure). So whatever q may be, the particle can go smoothly once it chooses 
the shortcut. There is no difference whether it returns back to segment 2 or not. 

In figure 6(c) and (d), when p = 0, the differences seem to be obvious, when we 
contrast these two figures, we are able to find that the extent of difference is decided 
by parameter q. While q changes from 0.1 to 0.9, the error changes from within 
6 X 10~^ to 0.015. This is simply because the smaller q means less particles to choose 
the shortcut. But if we compare figure 6(b) with 6(d), it can apparently be found that 
the error changes greatly when p is switched from 1 to 0, providing q is comparatively 
large. So we can say that it really matters whether a particle can return back or not 
when many particles choose the shortcut. If the particle can not return, a great traffic 
jam will happen at the certain site, refiected in figure 6(d) as a high-positioned blue 
star on the top of the figure. 

B. Simulations of advanced model 1 

Now we use numerical methods to simulate advanced model 1. In order to find 
how the phases of five segments change when parameters qi and q2 change, we set 
Pi = P2 = ^, thus corresponding to our equation analysis in section 3. We will see 
whether the numerical results correspond to our analytical results. 

Prom figure 7(a), we can see that all of the five segments are in low-density phase. 
The situation in segment 2 has nothing to do with the situation in segment 4. We 
can treat them as a connection of two basic models. The density rate in segment 
2 decreases when qi increases. The density rate in segment 4 decreases when q2 
increases. 

Prom figure 7(b), we can see that all of the five segments are in high-density phase. 
The situation in segment 2 has nothing to do with the situation in segment 4. We 
can treat them as a connection of two basic models. The density rate in segment 2 
increases when qi increases. The density rate in segment 4 increases when q2 increases. 

Prom figure 7(c) and 7(d), we can see that the situation in segment 2 and 4 are 
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FIG. 7: numerical simulations of advanced model 1 corresponding to different phases: (a)a = 0.3, (3 = 
0.8; (b)a = 0.8, 13 — 0.3; (c)a = /3 = 0.8, initial densities of all sites are high; (d)a = /3 = 0.8, initial 
densities of all sites are low. 
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FIG. 8: numerical simulations of advanced model 2 with insertion rate a = 0.3 and /3 = 0.8: (a)(72 
changes when qi remains constant; (b) qi changes when q2 remains constant. 

relevant to the initial density of each site. If the average initial density of all sites is 
high, then segment 2 and 4 will be both in high-density phase. If the average initial 
density of all sites is low, then segment 2 and 4 will be both in low-density phase. 



C. Simulation of advanced model 2 

Now we use numerical methods to simulate advanced model 2. In order to find 
how the phases of four segments change when parameters qi and q2 change, we set 
Pi = P2 = 1; thus corresponding to our equation analysis in chapter 3. We will 
see whether the numerical results correspond to our analytical results. Also, we are 
interested in whether there are interactions between the two parameters qi and q2- 
First, from figure 8 to figure 11, we see that the numerical simulation results 
are corresponding to our analytical results. The four segments should be all in low- 
density phase or high- density phase, and when segment 1 and 4 are in maximum 
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FIG. 9: numerical simulations of advanced model 2 with insertion rate a = 0.8 and /3 = 0.3: (a)(72 
changes when qi remains constant; (b) qi changes when q2 remains constant. 

current phase, the phase of segment 2 and 3 will be based upon the initial density of 
each lattice. 

If we focus on figure 8, we will find some interesting things by comparing figure 
8(a) and figure 8(b). In figure 8(a), we let q2 change but qi remain constant. We 
see that both the densities in segment 2 and segment 3 decrease in the same pace 
when 52 increases. In figure 8 (b), we let qi change when q2 remains constant. The 
result is that although the density in segment 2 decreases with increasing of gi, the 
density in segment 3 does not have apparent change. This phenomenon means that 
the change of qi can not influence the density in segment 3 greatly, but the change of 
q2 can really influence the density in segment 2 greatly. From this we know that the 
parameter q2 has priority to control the density of segment 2 and 3. 

From figure 9 to 11, if we compare (a) and (b) of each figure carefully, it is easy 
for us to find that for each figure, there exists the same phenomenon as in figure 8. 
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FIG. 10: numerical simulations of advanced model 2 with insertion rate a = 0.8 and /3 = 0.8 with 
high initial density for each lattice: (a)(3'2 changes when q\ remains constant; (b) q\ changes when 
^2 remains constant. 

In conclusion: in advanced model 2, the density of segment 2 and 3 depends on 
parameter 52 greatly. The change of will bring about the change of density of both 
of the segments. Once ^2 is set constant, parameter q\ will also has effects on the 
density in segment 2, but the density in segment 3 will remain in a relatively stable 
condition. 

V. CONCLUSION 

In this paper, TASEP with one or two shortcuts have been analyzed in details. 
We have found that when a particle arrives at the beginning of the shortcut, it faces 
the choice of whether to jump through the shortcut or not. If it chooses the shortcut, 
then a problem exists that whether it can return back to the ordinary road when it 
finds the road ahead blocked, according to which we have two different situations. 
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FIG. 11: numerical simulations of advanced model 2 with insertion rate a = 0.8 and /3 = 0.8 with 
low initial density for each lattice: (a)(?2 changes when q\ remains constant; (b) q\ changes when 
remains constant. 

After the study we know that if the particle can not return back, the beginning of the 
shortcut will most probably be filled with particles, producing a heavy traffic jam. 
This research offers us an idea that when we construct a road with a shortcut, the 
entrance of the shortcut, that is to say the place where the main road bifurcates, 
ought to be built widely. Once a driver chooses to jump the shortcut but find the 
signal showing "blocked ahead", he can immediately turn back to the ordinary road 
to prevent time-wasting. 

From the simulation of basic model, we find that no matter which situation is, the 
three segments should be in the following four phases: 

{L,L,L),{H,H,H),{M,L,M),{M, H, M) 

Whether they are in (M, L, M) or (M , H, M) is based on the initial density of the 
lattices. 
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Based on basic model, two advanced models have been studied. One of them 
have two shortcuts at different places while the other have two shortcuts at the same 
beginning. 

For advanced model 1 we have found that the five segments of the model should 
be in the following four phases: 

{L,L,L,L,L),{H,H,H,H,H),{M, H, M,H,M),{M,L, M, L, M) 

For advanced model 2 we have found that the four segments of the model should be 
in the following four phases: 

(L, L, L, L), (i/, H, H, H), (M, L, L, M), (M, H, H, M) 

From the simulation of advanced model 1, we have found that this model can also 
be regarded as two basic models, the segment 3 of the first being the segment 1 of 
the second. The probabilities of the choices of a particle facing the two shortcuts 
do not interfere with each other, which means that the probability of the choice of 
the particle facing one of the shortcuts does not infiuence the density in the other 
segment. This phenomenon leads us to suppose that if there are many shortcuts at 
different places of the road, we can view them as a connection of basic models. 

From the simulation of advanced model 2, we have found that the second shortcut 
can decide the density of both of the two segments in the middle while the first 
shortcut can only decide the density in segment 2. It offers us an idea that if we 
induce more drivers to go directly through the second shortcut, the density of the 
ordinary road will essentially decrease, providing the whole road is not so crowded. 

However, in this paper all of the proofs are under the condition p — 1. Future 
work can be focused on how to proof the phase situations for all p G [0, 1]. Also, two 
shortcuts overlapping each other will be studied in the future. 
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